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The recent WMAP data have confirmed that exotic dark matter together with the vacuum energy 
(cosmological constant) dominate in the flat Universe. Modern particle theories naturally provide 
viable cold dark matter candidates with masses in the GeV-TeV region. Supersymmetry provides the 
lightest supersymmetric particle (LSP), theories in extra dimensions supply the lightest Kaluza- Klein 
particle (LKP) etc. The nature of dark matter can only be unraveled only by its direct detection in 
q ■ the laboratory. All such candidates will be called WIMPs (Weakly Interacting Massive Particles) . In 

any case the direct dark matter search, which amounts to detecting the recoiling nucleus, following 
its collision with WIMP, is central to particle physics and cosmology. In this work we briefly review 
the theoretical elements relevant to the direct dark matter detection experiments, paying particular 
' attention to directional experiments, i.e experiments in which, not only the energy but the direction 

of the recoiling nucleus is observed. Since the direction of observation is fixed with respect the the 
earth, while the Earth is rotating around its axes, the directional experiment will sample different 
parts of the sky at different times during the day. So, since the event rates are different when looking 
at different parts of the sky, the observed signal in such experiments will exhibit very interesting 
Q^,) , and characteristic periodic diurnal variation. 
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PACS numbers: 95.35.+d, 12.60.J\ 



The combined MAXIMA-1 0, BOOMERANG 0, DASI and COBE/DMR Cosmic Microwave 
Background (CMB) observations imply that the Universe is flat and that most of the matter 
in the Universe is Dark 0. i.e. exotic. These results have been confirmed and improved by the 
recent WMAP data 0. Combining the the data of these quite precise experiments one finds: 

Q b = 0.0456 ±0.0015, Qcdm = 0.228 ± 0.013, ft A = 0.726 ± 0.015. 



Since the non exotic component cannot exceed 40% of the CDM [8| , exotic (non baryonic) matter is 
required and there is room for cold dark matter candidates or WIMPs ( Weakly Interacting Massive 
' Particles). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the 
observed rotational curves, see e.g the review 0, it is essential to directly detect such matter. Until 
dark matter is actually detected, we shall not be able to exclude the possibility that the rotation 
curves result from a modification of the laws of nature as we currently view them. This makes it 
imperative that we invest a maximum effort in attempting to detect dark matter whenever it is 
possible. Furthermore such a direct detection will also unravel the nature of the constituents of 
dark matter. The possibility of such detection, however, depends on the nature of the dark matter 
constituents. 

Supersymmetry naturally provides candidates for the dark matter constituents. In the most 
favored scenario of supersymmetry the LSP (Lightest Supersymmetric Particle) can be simply de- 
scribed as a Ma.jorana fermion, a linear combination of the neutral components of the gauginos and 
higgsinos [ill H2I IT3I [Til Ha . [lt| . In most calculations the neutralino is assumed to be primarily 



a gaugino, usually a bino. 

Since the WIMP's are expected to be very massive, mwiMP > 10 GeV, and extremely non 
relativistic, with average kinetic energy (T) w 50 KeV(mwiMp/100 GeV), they are not likely to 
excite the nucleus. So they can be directly detected mainly via the recoiling of a nucleus (A,Z) in 
elastic scattering. The event rate for such a process can be computed from the following ingredients: 
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An effective Lagrangian at the elementary particle (quark) level obtained in the framework 
of the prevai ling pa rticle theor y. For supersymmetry this is achieved as described, e.g., in 

Refs. uj& m Ea El El EI 0, 113, El- 



A well defined procedure for transforming the amplitude obtained using the previous effective 
Lag rangian from the quark to the nucleon level, i.e. a quark model for the nucleon, see e.g. 



iaEaMliS- This step n SUSY models is not trivial, since the obtained results depend 

crucially on the content of the nucleon in quarks other than u and d. 

Knowledge of the relevant nuclear matrix elements 2^, 22l , obtained with as reliable as possible 



many body nuclear wave functions. Fortunately in the case of the scalar coupling, which is 
viewed as the most important, the situation is a bit simpler, since then one needs only the 
nuclear form factor. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. 

We do not know for sure what this distribution is, but it is not expected to depend on the nature 
of the WIMP. In other words all WIMPs are expected to have the same velocity distribution 
and matter density. The particle density, however, which enters the event rate, is expected to 
be inversely proportional to the WIMP mass. In the present work we will consider a Maxwellian 
distribution. 

In the standard nuclear recoil experiments one has the problem that the reaction of interest does 
not have a characteristic feature to distinguish it from the background. So for the expected low 
counting rates the background is a formidable problem. Some special features of the LSP-nuclear 
interaction can be exploited to reduce the background problems. Such are: 

• The modulation effect. 

This yields a periodic signal due to the motion of the earth around the sun. Unfortunately 
this effect is small, < 2% for most targets. Furthermore one cannot exclude backgrounds with 
a seasonal variation. 

• Transitions to excited states. 

In this case one need not measure nuclear recoils, but the de-excitation 7 rays. This can 
happen only in vary special cases since the average WIMP energy is too low to excite the 
nucleus. It has, however, been found that in the special case of the target 127 I such a process 
is feasible [24| with branching ratios around 5%. 

• Detection of electrons produced during the WIMP-nucleus collision. 

It turns out, however, that this production peaks at very low energies. So only gaseous TPC 
detectors can reach the desired level of 100 eV. In such a case the number of electrons detected 
may exceed the number of recoils for a target with high Z 
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• Detection of hard X-rays produced when the inner shell holes are filled. 
It has been found 2jj that in the previous mechanism inner shell electrons can be ejected. 
These holes can be filled by the Auger process or X-ray emission. For a target like Xe these 
X-rays are in the 30 keV region with the rate of about 0.1 per recoil for a WIMP mass of 100 
GeV. 

In the present paper we will focus on the characteristic signatures of the WIMP nucleus interaction, 
which will manifest themselves in directional recoil experiments, i.e. experiments in which the 

direction of the recoiling nucleus is observed [H M H Im S3, H, 0, M \m ■ 

We will concentrate 

on the standard Maxwell-Boltzmann (M-B) distribution for the WIMPs of our galaxy and we will 



not be concerned with other distributions [221 [23, [2a, [40j, even though some of them may yield 



stronger directional signals. A mong those one should mention the late infall of dark matter into the 
galaxy, i.e caustic rings 41, 42l. I43l. l44l.l45t. dark matter orbiting the Sun 32 1 and Sagittarius dark 
matter [H]. 

We will present our results in such a fashion that they do not depend on the specific properties of 
the dark matter candidate, except that the candidate is cold and massive, tuwimp > 10 GeV. So, 
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since the nucleon cross section will be treated as a free parameter, the only parameters which count 
is the reduced mass, the nuclear form factor and the velocity distribution. So our results apply to all 
WIMPs. In a previo us pape r we have found that the observed rate is correlated with the direction 



of the sun's motion [34|, [35J, [36| . On top of this one will observe a relatively large time dependent 
variation of the rate due to the motion of the earth. Those features cannot be masked by any known 
background. To fully exploit these features the detector should be able to distinguish head from 
tail of the track, i.e. the sense of direction of motion 47[. An analysis of the requirements imposed 
on the directional detectors, focusing especially on energy threshold effects and minimization of 
background, has recently appeared [481 ] - The currently planned experiments, however, may not be 
able to distinguish the two possible senses along the line of nuclear recoil [49(. Such experiments 
cannot, e.g., measure the backward-forward asymmetry. 

Previous calculations considered the possibility of the variation of the event rate as the direction 
of observation of the nuclear recoil is varied in the galactic coordinates. This expected asymmetry 
is due to the motion of the sun with respect to the center of the galaxy. The apparatus will, of 
course, be oriented in a direction defined in the local frame, i.e. the line of observation will point to 
a point in the sky specified, in the equatorial system, by right ascension a and inclination 5. Thus 
due to the rotation of the earth the directional experiment will probe different parts of the galaxy 
during the day. Hence we expect to observe a diurnal variation. In light of this possibility we will 



utilise our previous directional calculations [34 , l35l . 
the previous results persist, as a diurnal variation. 



361 ] to explore, which characteristics, if any, of 



WIMP-NUCLEON CROSS SECTION 



The WIMP nucleon differential cross section in the non relativistic limit is given by 



da 



2?7K 



2m v 



2m r , 



(1) 



where p, p are the momenta of the initial, final WIMP and q the momentum transferred to the 
nucleon assumed to be initially at rest. Momentum and energy conservation are incorporated by 
suitable 8 functions. m x and m p are the masses of the WIMP and the nucleon respectively. M. n is 
the invariant amplitude at the nucleon level. 
Integrating over the final WIMP momentum we get 
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where v is the WIMP velocity, fj, p w m p is the reduced mass of the WIMP nucleon system and 
£ = p.q. The integration over q is trivial due to the (5 function. The allowed value of q is q = 2/x p u£. 
We thus get 



11,.. l2 (2^ p u£) At _ 1 | (1 i :2< 



da = - — \M n 

v Zir 



-d£ = —\M n \'4i4td£. 



u£ 2?r 1 

The scattering is forward, < £ < 1. We thus find the nucleon cross section to be 

1 



2tt 



\M n \ J 2i4 



(3) 



(4) 



In the present calculation we are not going to be concerned with calculating the nucleon cross section, 
which is perhaps the most important ingredient as was discussed in the introduction 

[13, EH EE 

14, 15, 3, U, HI, 3 2(3, 21 1 . We will instead treat it as an input parameter (see below). 
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WIMP-NUCLEUS CROSS SECTION 

Proceeding as above we obtain 

da = ±\M A \ 2 2fad£, (5) 

where now Ma is the corresponding amplitude at the nuclear level and u r is the WIMP-nucleus 
reduced mass. 

In this case instead of the variable £ it is advantageous to use the energy transfer Q = q 2 /(2mA)- 
We prefer to use the dimensionless variable u = (l/2)(qb) 2 , where b is the harmonic oscillator size 
parameter of the shell model of the nucleus. Thus 

u= l -( q bY = 2£v 2 e^i = ^^-- (6) 
£ V 2mu r vb 

Using the standard formula for the size parameter we can express u in terms of the energy transfer 
Q: 

u = 7T > Qo = (^Arrip)- 1 = 40A" 4 / 3 MeV , (7) 
Qo 



we thus get 



To proceed further we should relate |A^a| 2 to the more elementary amplitude |A^ n | 2 - This of, 
course, depends on the model. In the usual scace of the coherent scattering we get 

\M A \ 2 = A 2 F 2 (u)\M n \ 2 , (9) 
where F(u) is the nuclear form factor. We can thus relate the cross sections 

STANDARD (NON DIRECTIONAL) RATES 

The event rate is for a given WIMP velocity is given by 

cm = ^^A 2 F 2 ( u y n (!±) 2 V *L^, (ii) 

m x Am p VMp/ 2(p ir bv)' L 

where m t /Am p is the number of nuclei in a target of mass m t . Also 



d_R 

du m 



^-^V^AV n ( * V (12) 
m x Am p \ Hp J du 



where 

dt(u,v) _v 2 1_ 

du ~ V^ 2 ^ {U) 2{u r bv) 2 - [U) 

We must now fold dt(u,v)/du with the velocity distribution f(y,£) in the local frame. We will 
assume that the velocity distribution is Maxwell-Boltzmann in the galactic frame, namely 

f(y') = J-e-(y')\ y' = ^, (14) 

ITy/TT V 
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where v' is the WIMP velocity in the galactic frame and vq is the sun's velocity with respect to the 
center of the galaxy. We have, however, to transform the distribution it to the local frame: 

y' = y + z + S (sin ax — cos a cos jy + cos a sin 72) , (15) 

where 5 = vi/vq = 0.136 with V\ the velocity of the earth around the sun, a the phase of the earth 
(a=0 on June 3nd) and 7 ss 7r/6. 



No modulation (5 — 0) 

With S = we find the distribution 

ffat) = J- e -<v'+»ut+i). (16) 



Thus we get for the differential (with respect to the energy transfer) rate 

^ = ^ t>y ^ ^ = yi^^^, (it) 

au m x Am p \ fi p ) au au V 3 



^o(aV^) = [""' ydy 2n [° d£f(y, £), (18) 

where £ is a function of y and may come from the requirement that the WIMP velocity must not 
exceed the escape velocity. In practice £0 « 1 and also 

a = (V^firbvoy 1 . (19) 

The limits in y are obtained as follows: 

Vminiu) <^ £ = 1 =>■ y m in{u) = dy/u , y m(M ; = J/esc = 2.84, (20) 

where the escape velocity is t> esc = 2.84-u - 

We note that the function ^S>oi x ) depends only on the velocity distribution and is independent of 
nuclear physics (for £ — 1 wu l be- given below, Eq. (j2"B"|) ). It is a deceasing function of u as shown 
in Fig. [1] An additional suppression as the energy transfer increases comes, of course, from the 
nuclear form factor F(u). Note that the nuclear dependence of the differential rate comes not only 
from the form factor, but via the parameter a as well. Integrating the differential event rate from 
Umin = Eth/Qo, which depends on the energy threshold, to 

Umax — Uesc 

I a 2 we obtain the total rate 



n-^m'A% r-d.* <*»> 



m x Am p VMp/ Ju in 



The modulation effect 

Now the velocity distribution in the local frame becomes more complicated. We get 

f(y, £, d>, a, 7, 6) = J- e -(v 2 +^+b+c) (22 ) 

with 

A = 5 cos a sin 7 + 1, B = 5 2 + 2 cos a sin j5 + 1 , (23) 



G 




and 

C = 2yS yl — £ 2 (sin a cos <fi — cos a cos 7 sin <j>) . (24) 

We now notice that to leading order in S the term C does not contribute, since the integration over 
(f> yields zero. Then the integration over <fi and £ can be done analytically to yield 

2e-« 2 - 52 - 2l5cosQsin ^- 1 sinh ( 2 (8 cos a sin jy + y)) 

f(y, a, 7, S) = rTk • i — ^ • ( 25 ) 

y 7r(d cos a sin jy + y) 

The integration over y: 

_ rVe.sc 

^(a-y/u, a, 7, (5, y esc ) = / ydyf(y, a, 7, 5), 
can also be done analytically to yield 

g — 5 2 — 2 cos a sin "y<5+(<5 cos a sin 7 + I) 2 — 1 

*(x,a,7,(5,?/ esc ) = — : — 

2{o cos asm 7 + 1) 

erf(— x + S cos a sin 7 + 1) + erf (a; + S cos a sin 7 + 1) 
— erf (5 cos a sin 7 — y csc + 1) — erf (<5 cos a sin 7 + y csc + 1) , (26) 

with x = a^/u and erf(t) the well known error function. Expanding this in a power series up to 
linear terms in 5 we get: 

= * (x) + *i(x,7, <5) cos a, (27) 
*o(a:) = 7, 0, y esc ) = ^ ( crf(l - x) + erf(l + z) - erf(l - j/ esc ) - erf(l + y esc )) , (28) 
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FIG. 2: The function ^\{x, 1/2,0.135) arising from the velocity distribution, which enters the 
modulated differential rate. Notice the change in sign. The position this happens depends on the 
reduced WIMP-nucleus mass. This dictates whether the modulation of the total rate is positive or 

negative. 



for the zeroth order term and 



1 



*i(a,7,<y) = -5 



erf(l - a:) - erf(a; + 1) + erf (1 - y c sc ) + erf (y csc + 1) + 

2e -(x+l) 2 2e - fecsc -l)2 2e _ (yoBc + i) 



2e-(*- 1 ) 2 



sin(7), 



(29) 



for the coefficient of cos a of the first order term. This is exhibited in Fig. [5] We thus proceed as 
above with 



(30) 



Thus we finally get 



dt dr _ dt dh 
du du du du 



t — » t(l + /icosa), h 



du- 



dh 
du 



R=£2L^[^L \ AU(l + hcosa) 
m x Am p \ fi p ' 



(31) 
(32) 

(33) 



DIRECTIONAL RATES 

In this instance the experiments will attempt to measure not only the energy, but the direction of 
the recoiling nucleus as well [36| , • Let us indicate the direction observation by: 



e = (e X7 e y ,e z ) — (sin 9 cos $, sin 9 sin <E>, cos 9). 



(34) 
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Now (|13l) must be notified since the variables u and v are no longer independent. Furthermore no 
integration over the azimuthal angle specifying the direction of the outgoing nucleus is performed. 
Eqs. (fT2")) and ((13] ) now become 



dR 

du 
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Px m t 
m x Am p 



f dt(u, v) 
\ du 



V< v 2 >A 2 a n 



F 2 (u) 



1 fdt{u,v) 
2^ V du 



dir 
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V< V 2 > 



2{p r bv) 2 



v/v 



(35) 



(36) 



The factor 1 / (2tt) enters since we are going to make use of the same nuclear cross section as in the 
non directional case. 

Because of this constraint and the fact that the limits of integration over the WIMP velocity depend 
on the angles and We thus find it convenient to use a more convenient system (X,Y,Z), in 
which the polar axes Z is in the direction of the recoiling nucleus. The needed transformation is 









-I 








(37) 



(38) 



cos(0) cos($) - sin($) cos($) sin(0) 
= ( cos(0) sin($) cos($) sin(0) sin($) 
- sin(0) cos(0) 

and thus the three components of the velocity are given by: 

vx\ /cos(0)cos(3>) cos(0)sin($) - sin(0) \ / <5sin(a) 
vy — sin($) cos($) — 8 cos(a) cos(7) 

vz ) \ cos($) sin(0) sin(0) sin(<f>) cos(0) J \ £cos(a) sin (7) + 1 

8 cos(0) cos($) sin(a) — (5 cos(a) sin(7) + 1) sin(0) — 8 cos(a) cos(7) cos(0) sin(<I>) 
—6 cos(a) cos(7) cos(<£>) — 8 sin(a) sin(<I>) 

cos(0)((5cos(a) sin(7) + 1) + (5cos($) sin(a) sin(0) — 8 cos(a) cos(7) sin(0) sin($) 
The exponential of the M-B distribution ( Eq. (TT4"|) ) in this new local frame becomes 

y 2 + 2y - 5y/i - £ 2 sin(^) ( cos(a) cos(7) cos($) + sin(a) sin($)^ + £ ( cos(0)(<5 cos(a) sin(7) + 1) 

+ <5sin(0)(cos(<I>) sin(a) — cos(a) cos(7) sin($)) J — \J\ — £ 2 cos(^>) ( (8 cos(a) sin (7) + l)sin(0) 

+ 8 cos(0)(cos(a) cos(7) sin($) — cos ($) sin(a)) \ J + 8 2 + 28cos(a) sin( 7 ) + 1, (39) 

where in this new frame the direction of the velocity is given by (9, </>) and v.e — cos 9 = £. 
Expanding in powers of 8 up to first order we get 



f(y, e, 0, 0, $, a, 7, 8) = e -(a 2 +2^cos(e)-2,yr=^cos W )si„(e)+i) ( 1+pc0+ps0+pcc+pcs+psc+pss ^ 

(40) 

with 



pcO = 8 



pec = 8 



2 sin(7) — 2y£ ( cos(0) sin(7) — cos(7) sin(0) sin($) 
psO — 8 ( — 2y^cos($) sin(0)^j sin a, 
2yy/l — £ 2 ( sin(7) sin(0) + cos(7) cos(0) sin($) N 
pes — 8 \ 2yy/l — £ 2 cos(7) cos(<£>)^ cos a sin 0, 



cos a, 



cos a cos < 



(41) 
(42) 
(43) 
(44) 



psc 
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S (—2y\/l - £ 2 cos(6) cos($)^ sin a cos 0, (45) 
pss = 5 \ 2yy/l — £ 2 sin($)^ sin a sin 0. (46) 



We now find that the integrals over <fi can be done analytically 

flit p2TT 



e acos4> d(j)= / e «"in*d0 = 27rJ (o), (47) 
'0 Jo 

<-27T /.2-ir 



Z7T />Z7T 

/ e QCOS<#, cos(/)d0= / e Qsin *sin0d< ? !) = 27r/ 1 (a), 
Jo Jo 



(48) 



2 IT 



e asin0 cos0d0 = O. (49) 



/o Jo 
where Io(a) and Ii(a) are the modified Bessel functions. 

No modulation (5 = 

After the <j> integration we get 

2 



f(y,Z,Q) = -^e-v'-^c-CB)-!^ ( 22/ yrr^ sin( e)) . (50) 
We can now use the 5 function to perform the integral over £ to get: 

f(x, y, 9) = JL e -y 2 -2-os(e)-i /o ( 2 V^^sin(e)) , (51) 

where x = a^/u. The function ^f r (a^u, 9, $) can only be obtained numerically: 

6) = |//(oV«, l/> ©)dtf. (52) 

J a^/U 

Once this function is known the computation of the event rate proceeds in a fashion analogous to 
the standard (non directional) case. The behavior of this function is shown in Fig. [3l 

Modulated Rates in Directional experiments 

It is clear that the terms pes and pss make no contribution. Thus after the (f> integration we get 
f(y, 9, $, a, 7 , 8) = -l e -v 2 -^»(e)-i {qcQ + qgQ + qcc + qss) (53) 

\/7T 



where 
qcO = S 



2 sin 



in(7) - 2y£ ( cos(0) sin(7) - cos(7) sin(6) sin($)^ I ^2yy/l - £ 2 sin(9)) cos a, 

(54) 



qs0 = 5[ - 2y£cos($)sin(9)) I [2yy/l-£ 2 sin(9)) sin a, (55) 



qcc = S 



Wl ~£ 2 ( sin( 7 ) sin(6) + cos( 7 ) cos(9) sin($)) h (tyy/l - £ 2 sin(9)) cos a, (56) 



qsc = 5 (-2yy/l - £ 2 cos(9) cos($)) /i V 1 ~ C 2 sin(9)) sin a. (57) 
Proceeding as above we can use the 6 function to integrate over £ and get two modulation amplitudes: 
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X — CLy/u — > 

FIG. 3: The function (x, 0) entering the directional differential rate. Here we exhibit the 
curves corresponding from zero to it in steps of 7r/10 ( in the order from down up). We see that, 
with the possible exception of very small x = a^fu (small energy transfers) the differential rate 
increases with the angle (maximum opposite to the sun's direction of motion). 



• The amplitude that multiplies the cos a term is: 



with 
rcO = 5 



Ux, y, 0, $, S) = — e-J/ 2 - 2 *^ )- 1 ^ + rcc), (58) 

2 sin(7) - 2x ( cos(0) sin(7) - cos(7) sin(0) sin(<I>)^ I Q (2\/ y 2 - x 2 sin(0)^ , 

(59) 

rcc = 2<5vV - x 2 ( sin( 7 ) sin(0) + cos( 7 ) cos(0) sin($)) I x (2^y 2 - x 2 sin(0)) . (60) 
Thus we numerically get 

¥* r (ay^, 0, $, S) = f"" y / c (oVS, y, 0, M)dy. (61) 
• The amplitude that multiplies the sin a term is of the form: 



f s (x, y, 0, <f, 6) = e V-to«(e)-i (r(0 + rsc); 

V7T 



where 



rsO = -28x ( cos($) sin(0) I ( 2 y/y 2 ^ 2 sin(0) 



Then finally: 



-2^Vy2 - x 2 ^ cos (0)cos($)) Ji ^v/y 2 -x 2 sin(6)) 



(62) 

(63) 
(64) 

(65) 
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X — dy/u — > 



FIG. 4: The function 9? r {x, 6, $, 5) on the left and 6, $, (5) entering the modulated 

directional differential rate as coefficients of the cos a and sin a respectively . They correspond to 
O = (3/10)7r, i.e. almost parallel to the sun's direction of motion. The overall rate in this direction 

is expected to be small. The values of $ range from zero to 2n in steps of it/ 5 ( the style of the 
curves is the same as in Fig. [31 but the step here is n/5). We see that the $ dependence is strong 
at small energy transfers and but it gets smaller as the energy transfer increases. 




FIG. 5: The same as in Fig. 2] for = ir/2, i.e in a plane perpendicular to the sun's direction of 

motion. 



The differential rate is given: 



dR\ p v m* I U r \ .o /2 



du Jdir m X Am P V/W V 3 

-!- f 9f r (ay/u, 6, $) + r (aV5, 6, $, S) cosa + *f r (a^, 6, $, <5) sina) . (66) 

Z7T V / 

The behavior of the functions *f r (x, 6, <5) and *f ""(oh/u, 6, $, <S), x = ay^, is exhibited in 
Figs H] -[61 The total event rate is obtained after integrating the above expression over the energy 
transfer. It can be cast in the form: 



p Px m t 



\ 2 1 

— J A 2 —tdir(@, $) ( 1 + h c (Q, $) cosa + h s (Q, $) sin a) . (67) 



We found it convenient to use a new relative parameter k (to avoid or minimize the dependence on 
the many parameters of the particle model) as well as the parameters h m and «o given by: 

k(8,$) = — , /i ot (6,$)cos(q; + q;o) = ft c (8 1 $)cosa + ft a (e,$)sina. (68) 
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FIG. 6: The same as in Fig. 0]for 9 = (9/10)7r, i.e very close to the direction opposite to that of 
the sun's velocity. In this case the two coefficients are almost identical and, for given $, the shape 
is similar to the non directional case (see Fig. [5]). 



These parameters depend, of course, on the reduced mass /x r , the WIMP velocity and to some 
extend on the nuclear physics via the nuclear form factor. The parameter n gives the retardation 
factor of the directional rate, over and above the factor of l/(2w), compared to the standard rate. 
Since, however, the unmodulated amplitude is independent of $ one can integrate over <E> so that 
the suppression of l/(27r) drops out for this term. Because of the existence of both cos a and sin a 
terms the time dependence will be of the form cos(a + ao), with the phase ckq being direction 
dependent. Thus the time of the maximum and minimum will depend on the direction. In other 
words the seasonal dependence will depend on the direction of observation. So it cannot be masked 
by irrelevant seasonal effects. 

Before concluding this section we like to consider the case connected with the partly directional 
experiments, i.e. experiments which can detemine the line along which the nucleus is recoiling, but 
not the sense of direction on it. The results in this case can be obtained by summing up the events 
in both directions, i.e. those specified by (9, $) as well as (ir — 9, $ + 71"). A given line of observation 
is now specified by 9, <I> in the range: 

< 9 < vr/2, 0<$<7r. 



SOME APPLICATIONS 



In this section we are going to apply the formalism of the previous section in the case of two 
popular targets: i) The light target 32 S appearing in CS2 involved in DRIFT (2^1 and ii) The heavy 
targets 127 I , which has been employed in the DAMA experiment 
by the XENGN collaboration. We will not consider energy thresh hold effects and we will ignore 
quenching factor effects. We will consider only the coherent mode, but the results obtained for the 
functions t, n are not expected to be radically modified, if one considers the spin mode (the quantities 
h m and a$ will be discussed in a future publication. The nuclear form factor employed was obtained 
in the shell model description of the target and is shown in Figs. [7] and [5] for a light and a heavy 
target respectively. 



51, 51 and 131 Xe 53, 54 



employed 



The non directional case 



The parameter t and h entering the non directional case is shown in Figs. [9] and 1101 Using the 
above parameters t and employing Eq. (|33[) we obtain the time independent total event rate. The 
obtained results are shown in Fig. [11] The total modulation, since it is defined relative to the time 
independent part, is still given by Fig. [101 
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Q — y keV 



(a) (b) 

FIG. 7: (a) The form factor F 2 {u) for 32 S employed in our calculation with u = Q/Qo, Q the 
energy transfer to the nucleus and Qq = 404 keV. (b) The same quantity a function of the energy 

transfer Q. 




(a) (b) 

FIG. 8: The same as in Fig. [7] in the case of a heavy target (A=127 or A=131). 




300 400 500 



(a) (b) 
m x — > GeV 



FIG. 9: The quantity t co h, as a function of the WIMP mass in GeV, is shown for Q m in — for a 

light system, 32 S (a) and m I or 131 Xe (b). 



The directional case 
The parameter k entering the directional case is shown in Figs. [T2]fT5l 

DIURNAL VARIATION 

Up to now we have considered the event rate in a directional experiment in fixed direction with 
respect to the galaxy in the galactic system discussed above. The apparatus, of course, will be 
oriented in a direction specified in the local frame, e.g. by a point in the sky specified, in the 
equatorial system, by right ascension a and inclination S. This will lead to a diurnal variation of 
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(a) 



j>oo ft? 



-400- _500 



GeV 



(b) 



FIG. 11: The coherent total event rate, as a function of the WIMP mass in GeV, is exhibited 
assuming a nucleon cross section a N = 10~ 32 cm 2 for Q m i n = for a light system, 32 S (a) and a 

heavy system, 131 Xe (b). 




(a) (b) 
Q — > radians 

FIG. 12: The parameter k (see text) in the case of the target 32 S for a WIMP mass of lGeV, in the 
case the sense of direction is known (a) and when the sense of direction cannot be distinguished (b) . 



the event rate [55j |. 

The galactic frame, in the so called J2000 system, is defined by the galactic pole with as- 
cension ax = 12 h 51 m 26.282 s and inclination Si = +27° 1 42. Ol" and the galactic center at 
a 2 = 17 h 45 m 37.224 s , S 2 = -(28° 56' 10.23"). Thus the galactic unit vector y, specified by 
(ai,5i), and the unit vector s, specified by (a2,^2), can be expressed in terms of the celestial unit 
vectors i (beginning of measuring the right ascension), j and k (the axis of the Earth's rotation). 
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FIG. 13: The same as in Fig. 



(b) 

radians 

2] for a WIMP mass of 100 GeV. 




0.25 0.5 X- 75 1 1 ■ 25 1 ■ 5 



S£ 0.7 - 



(b) 



G — ^ radians 

FIG. 14: The same as in Fig. [T21for a a heavy target. 
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(a) 



(b) 



O — > radians 

FIG. 15: The same as in Fig. [TJfor a WIMP mass of 100 GeV. 



One finds 

y = -0.868i - 0.198.5' + 0.456A; (galactic axis) , 

x = —s = 0.055i + 0.873j + 0.483fc (radially out towards the sun) , 

z = x x y = 0.494i - 0.445 j + 0.747A: (the sun's direction of motion). (69) 

Note in our system the x-axis is opposite to the s-axis used by the astronomers. Thus a vector 
oriented by (a, 5) in the laboratory is given in the galactic frame by a unit vector with components: 

y \ / —0.868 cos a cos S — 0.198 sina cos <5 + 0.456 sin<5 \ 
x = 0.055 cos a cos S + 0.873 sin a cos S + 0.4831 sin S . (70) 
zj \ 0.494 cos a cos 5 -0.445 sin a cos 8 + 0.747 sin S J 
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This can also be parameterized [56[ as: 

S7C0S<5 cos (a — ao) — sin7 ( cosS cos (Op) sin (a — olq) + sin 8 sin (0p)) 



.r = cos ' 



y = cos (Op) sin (5 — cos <5 sin (a — ao) sin (Op) , 
z = cos 5 cos (a — ao) sin 7 + cos 7 ( cos 5 cos (0p) sin (a — ao) + sim5 sin (#p)) , 



(71) 
(72) 
(73) 



where ao = 282.25° is the right ascension of the equinox, 7 w 33° was given above and #p = 62.6° 
is the angle the Earth's north pole forms with the axis of the galaxy. Due to the Earth's rotation 
this unit vector, with a suitable choice of the initial time, is changing as a function of time 

x = cos 7 cos 5 cos (^~7jT^j ~ sin 7 f cos 8 cos (Op) sin + sm ^ sm (0p)\ > (74) 



2nt 

y = cos (Op) sin <5 — cos 8 sin ( ) sin (Op) , 



(75) 



cos 1 -jr^j c °s (5 sin 7 + cos 7 ^cos 8 cos (6*p) sin ) f ^iiui sin ( W;.. ) 



(76) 



where T is the period of the Earth's rotation. In fact the angles 8 and $ discussed above are given 
by 



6 = cos 1 . 



(77) 



$ = If 



> 0,lf 



y > 0, tan 



, tan 



3tt 
T 



.If 



y > 0, tan 



— , tan 

2' 



(78) 



$ = H( 



tan 



r(*)ff(y) (tan- 1 (HI)) +H(x)H(-y) 
+ H(-x)H(y) (tan- 1 (| ||) + |) + H(-x)H(-y) (tan" 1 
Some points of interest are: 



3_7T 

T 

y 

a; 



The celestial pole: (y,x,z) = (0.460, 0.484, 0.745) => 
The ecliptic pole: (y,x,z) = (0.497,0.096,0.863)^- 
The equinox: (y,x,z) = (-0.868,0.055,0.494): 



= 62.6°, = 57°), 
= 62.6°, = 83.7°), 
(0 = 150.2°, <£= 83.7°). 



(79) 



(80) 



where is defined with respect to the polar axis (here y) and (p is measured from the x axis towards 
the z axes. 

The angle scanned by the direction of observation is shown, for various inclinations 8, in Fig. 16. 
We see that for negative inclinations, the angle O can take values near tt, i.e. opposite to the direction 
of the sun's velocity, where the rate attains its maximum. Additionally, the angle $ which exhibits 
rather complicated behavior, is shown for various inclinations 8 in Fig. 17. 

Some typical paths of the orbit of the orientation point in the galactic coordinate system is shown 
for various inclinations in Fig. 1181 

The equipment scans different parts of the galactic sky, i.e. observes different angles 8. So the rate 
will change with time depending on whether the sense of the recoiling nucleus can be determined 
along the line of recoil. The results depend, of course, on the WIMP mass and the target employed. 
We will consider a light and an intermediate-heavy target. 
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FIG. 16: Due to the diurnal motion of the Earth different angles 6 in galactic coordinates are 
sampled as the earth rotates. The angle 9 scanned by the direction of observation is shown for 
various inclinations 5. The intermediate thickness, the short dash, the long dash, the fine line, the 
long-short dash, the short-long-short dash and the thick line correspond to inclination 
5 = — 7r/2, — 37r/10, — 7t/10, 0, n/10, 37r/10 and n/2 respectively. We see that, for negative 
inclinations, the angle can take values near tt, i.e. opposite to the direction of the sun's velocity, 

where the rate attains its maximum. 




a — > 



FIG. 17: The same as in Fig. [16] in the case of the angle <&. 



The case of a light target 

The time dependence of K is exhibited in Fig. [19] 

The case of an intermediate-heavy target 

The time dependence of k is exhibited in Fig [23 
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FIG. 18: The circular path followed by the point of the direction of observation as seen in the 
galactic system due to the Earth's rotation for various inclinations 5. The galactic axis is indicated 

upward. For the path notation see Fig. [TBI 



DISCUSSION 



We have seen that, since the expected event rates are small, which means that the background 
problems become formidable, one should exploit all the signatures provided by the reaction. One 
such signature is provided by the modulation effect h. Unfortunately, however, this effect is small. 
Furthermore, for heavy targets, its sign, which determines the location of the maximum, is not 
certain. It depends on the unknown WIMP mass. 

Better signatures are expected in directional experiments in which one measures not only the 
recoil energy but the direction of the nuclear recoil as well. Some of the requirements that should 
be met by such detectors have recently been discussed 47, 48|. To fully exploit the advantages of 
such detectors one should be able to distinguish between recoils with momenta p and — p (sense of 
direction), which now appears to be feasible. Some of the predicted interesting features of directional 
event rates persist, even if it turns out that the sense of motion of recoils along their line of motion 
cannot be measured. 

Such experiments given a sufficient number of events provide an excellent signature to discriminate 
against background. One, of course, gets a smaller rate by observing in a given direction. In the 
most favored direction, opposite to the sun's direction of motion, the event rate is ss down from 
that of the standard non directional experiments, if a specific angle $ is chosen. Since, however, k is 
independent of the angle one can integrate over all azymouthal and thus the retardation is only 

KRJl. 

Finally we have seen that in directional experiments the relative event rate for detecting WIMPs 
within our galaxy, as given by the parameter k, will show a periodic diurnal variation due to the 
rotation of the Earth. The parameter k is essentially independent of any particle model parameters 
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t/T — > 

FIG. 19: The time dependence (in units of the Earth's rotation period) of the parameter k for 
various inclinations S in the case of the target CS2. Due to the diurnal motion of the Earth 
different angles O are sampled as the earth rotates. For the notation see Figs [T5l and [TBI 



other than the WIMP mass. It does depend on the assumed velocity distribution. The time variation 
is larger in the case of light WIMP and/or light target. In all cases it is preferable to select a direction 
of observation with negative inclination, since, then, both the absolute value of k and its variation 
are larger. The time variation is larger in fully directional experiments. But, even if the sense along 
the line of motion is not observed, one will observe a time variation, but it will be smaller. This is 
hardly surprising, however, since, then, the observed event rate is the average of the two senses of 
direction. Note that in this case the time required to come to the same point is half of that of the 
previous case (see Figs [T9ll20|) . 

The time variation arising after the inclusion of the modulation parameters (h m ,a) or (h c ,h s ) 
is expected to be even more complicated, since these parameters depend on both O and <J>. One 
expects a diurnal variation on top of the annual variation characteristic of the usual modulation effect 
entering both directional and non directional experiments. Such effects will be discussed elsewhere. 
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t/T — > 

FIG. 20: The same as in Fig. [TO] in the case of the Iodine target. 
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